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Utilizing recent neutrino dimuon production measurement from NuTeV the assumptions on 
the determination of the strangeness content of the nucleon within the dynamical approach 
to parton distributions are investigated. The data are found to be in good agreement with 
the predictions derived from our (GJR08) dynamical parton distributions, which have been 
generated entirely radiatively starting from vanishing strange input distributions at an optimally 
chosen low resolution scale. Further, the data induce an asymmetry in the strange sea which is 
found to be small and positive in agreement with previous results. 

Parton distributions and their implications have been recently studied within the dynamical parton model 
approach [1,2]. Since the dynamical parton distributions at Q 2 > 1 GeV 2 are QCD radiatively generated 
from valencelike 1 positive definite input distributions at an optimally determined low input scale Qq < 1 
GeV 2 , the predicted steep small Bjorken-x behavior of structure functions is mainly due to QCD-dynamics 
at x < 10~ 2 . Alternatively, in the common "standard" approach the input scale is fixed at some arbitrarily 
chosen Qq> 1 GeV 2 and the corresponding input distributions are less restricted, for example, the observed 
steep small-x behavior of structure functions and consequently of the gluon and sea distributions has to 
be fitted here. Furthermore the associated uncertainties encountered in the determination of the parton 
distributions turn out to be larger, particularly in the small-x region, than in the more restricted dynamical 
radiative approach where, moreover, the evolution distance (starting at go < 1 GeV 2 ) is sizably larger (see 
[1,2] and references therein). 

As in previous dynamical determinations [3, 4], since the data sets used are insensitive to the specific 
choice of the strange quark distributions, the strange densities of the dynamical distributions in [1, 2] have 
been generated entirely radiatively starting from vanishing strange input distributions: 

s(x,Q 2 )=s(x,Q 2 )=0 (1) 

at the low input scale 2 . In order to investigate the plausibility of these assumptions, we confront here 
predictions derived from dynamical distributions determined in this way with data which are particularly 
sensitive to the strangeness content of the nucleon. 

For this purpose we have chosen the latest and most precise measurements of neutrino dimuon pro- 
duction from V^- and v^-iron deep inelastic scattering (DIS) interactions of NuTeV [5]. The measured 

'Valencelike refers to ay > for all input distributions xf(x, Qq) « x"f(l — x) b f, i.e., not only the valence but also the sea and 
gluon input densities vanish at small x. 

2 In the "standard" case, where 2o > 1 GeV 2 , the strange input distributions were chosen s(x,Qq) = s(x, Qq) = \{u{x, Qq) + 
d(x,Q^)), as is conventional [1, 2]. 



cross-section requires the muon from the semileptonic charm decay to have an energy greater than 5 GeV, 
therefore the theoretical predictions need to be corrected from detector acceptance; the "forward" cross- 
section is given by [5]: 

— (x,y,E m ) = B c £f(x,y,E m ) -j-j- {x,y,E v(9) ) , (2) 

where srf denotes the acceptance corrections, and d< f x ^y is the inclusive DIS cross-section for charged 
current charm production; fragmentation effects are included in the acceptance corrections as have been 
calculated in [5, 6] at NLO of QCD. The cross-section for (anti)neutrino charm production has been given 
up to NLO in [7, 8] within the "fixed flavor number scheme" (FFNS), i.e., besides the gluon, only the light 
quark flavors (u, d,s) should be included as (massless) partons within the nucleon (for a recent discussion 
on the treatment of heavy-quark masses see, for instance, [9]); for consistency we will therefore use our 
(GJR08) NLO FFNS dynamical parton distributions of [1] for the present investigations 3 . 

Since the NuTeV collaboration used an iron target, the mass M in Eq. (2) denotes the average nucleon 
mass and, moreover, the parton distributions have to be corrected for nuclear effects; we use the NLO 
corrections of [10]. These nuclear corrections were obtained using the previous set of NLO FFNS dy- 
namical parton distributions of GRV98 [4] (which are very similar to our GJR08, see [1]) and therefore 
are especially suited for combination with our GJR08 NLO FFNS. As is clear from Eq. (2), the absolute 
normalization of the measurements is directly related to the semileptonic branching ratio B c , hence we will 
allow its value to float within the experimental error of the value used by the experiment B c = 0.099 ±0.012 
[5], as is common for the (fully correlated) normalization errors in global QCD analyses. However, all our 
distributions prefer the highest allowed value. 

The results obtained with the GJR08 distributions are shown in Fig.l, where a good agreement is 
observed; the % 2 obtained is 65 for 90 data points 4 . This agreement demonstrates the compatibility of the 
data with the conditions of Eq. (1) and shows that in the dynamical case, where the input distributions are 
parametrized at an optimally chosen low input scale go =0.5 GeV 2 [1], the strange sea can be generated 
entirely radiatively starting from: 

s + (x,Ql)=s(x,Q 2 Q )+s( X ,Q 2 )=0 (3) 

In addition, the small differences between neutrino and antineutrino data are well known to induce a 
small asymmetry in the strange sea [5, 12-17]. In order to evaluate this asymmetry within our framework, 
we parametrize a new input distribution: 

s-(jc,fig) = s{x,Q 2 )-s(x,Q 2 )=Nx a (\ -x) b {\ - |) (4) 

where Qq = 0.5 GeV 2 is (fixed to) the input scale of our GJR08 NLO fit [1] and the function is constrained 
by the quark-number sum rule /J dx s~{x, Q^) — 0. Eqs. (3)-(4) imply that the strange input distributions, 
s(x,Qq) and s(x,Qq), will in turn be negative (positive) at the input scale, by construction. This is not 
a problem as long as at perturbative scales, say for Q 2 > 1 GeV 2 , both strange distributions s(x,Q 2 ) and 
s(x,Q 2 ) become manifestly positive due to the QCD evolution, as is the case. 

We have tried several (more elaborated) parametrizations without finding any improvement. As a matter 
of fact, even this "minimal" parametrization showing the required properties (vanishing at small and large 
x and crossing the x-axis at least once) is too flexible for the restricted x-range of the data that we are 
using (0.02 < x < 0.3). In practice there is a high correlation between a and b, so that both parameters 
should not be kept free simultaneously, in particular for the error estimation. We have chosen to fix b = 25, 
which is comparable to the value obtained for the analogous parameter of the light-quark sea-asymmetry 
input distribution bj_ u ~ 17 [1], but a bit larger as would be expected for the (heavier) strange quark; 
several values in the range 20 to 30 give similar results. This leaves us with 2 free parameters for the 

3 The use of the expressions in [8] together with "variable flavor number scheme" distributions, e.g. in [5], is inconsistent. 
4 Since the NuTeV collaboration has inflated the errors as a way of (not) taking into account the point to point correlations, this 
level of agreement briefly corresponds to 1 cr [11]. 
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Figure 1 : Comparison of the NuTeV "forward" dimuon cross-section data [5] with predictions using the (strange- 
symmetric) GJR08 distributions as well as our newly determined "asymmetric" results. 
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Figure 2: Out result for the strangeness asymmetry in the nucleon at Q 2 = 16 GeV 2 appropriate for the NuTeV 
experiment (cf. Fig. 3 of [5]). The results of MSTW2008 [16] are also shown for comparison. 

strangeness asymmetry; the crossing point of the input asymmetry xq is determined through the no-net- 
strangeness condition. The optimal values of these parameters have been determined following [1]; we get 
JV = -0.017±0.016 and a = 0.22 ±0.09. 

Since the rest of the data used in our global fit [1] are essentially insensitive to an asymmetry in the 
strange sea 5 and, moreover, we continue to generate the strange sea entirely radiatively through Eq. (3), the 
parameters quoted in [1] remain unchanged and our present result for the asymmetry can be used together 
with the GJR08 distributions. The evaluation of uncertainties is performed as in [1], however, since the 
free parameters introduced here are completely uncorrelated with the ones determined in [1], and there is 
no compatibility problems with the data, we do not find necessary to increase the tolerance parameter and 
use T = 1 for the 1 a uncertainties of s~ (x, Q 2 ). 

The newly obtained "asymmetric" results are also shown in Fig. 1, where it can be seen that the 
(anti)neutrino data prefer (smaller)larger values, i.e. the data favor an increase of the s distribution and 
a decrease of the s, in other words, a positive asymmetry in the relevant 0.01 < x < 0.1 region; for x val- 
ues larger than about 0.1 no significant changes are appreciated. After introducing the asymmetry, the 
X 2 value improves to 60 for 90 data points, although the predictions from the strange-symmetric GJR08 
distributions are rather similar and the differences lie within the 1 a bands. 

Our result for the strangeness asymmetry in the nucleon are shown in Fig. 2 at Q 2 — 16 GeV 2 appropriate 
for the NuTeV experiment, and can be directly compared with Fig. 3 of [5]. Although due to the large errors 
both results are in general agreement, the peak of our asymmetry is lower and placed at a slightly smaller 
value of x. The results of MSTW2008 [16] are also shown in Fig. 2 for comparison. They are rather 
similar in size to ours, despite the fact that in [16] older (and less precise) data from CCFR [11] have also 
been included and this tends to reduce the asymmetry [15]. Note, however, that our asymmetry is much 
more suppressed for large x> 0.2, where the data are in excellent agreement with our (strange-symmetric) 
GJR08 distributions [1], as can be seen in Fig. 1. 

'Actually, there are some marginal dependences, e.g. in or in the Drell-Yan cross-sections, which are completely negligible, 
for example, their effect in x 2 is well within the rounding error of the values quoted in [1]. 
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The changes in the strange-asymmetric distributions as compared with the original GJR08 are rather 
small, e.g. at g 2 = 100GeV 2 they reach at most 5% in the relevant lCT 3 <x<0.3 region, and are comparable 
with the uncertainties in the distributions, which are of a few percent as well. This being the case, the 
original strange-symmetric GJR08 distributions should suffice for most applications, moreover since most 
observables depend essentially only on s + (x, Q 2 ). Furthermore, since we continue to generate the strange 
distributions radiatively starting from Eq. (3), the increase in the uncertainties encountered in common 
"standard" fits, where s + (x, Qq) has to be fitted, is avoided in the more constrained dynamical case, which 
uncertainties should be very similar to the ones of GJR08 in most cases. 

The strangeness asymmetry is however relevant for applications especially sensitive to the strange 
content of the nucleon, as has been shown, for instance, in relation with the so-called the "NuTe V anomaly" 
(see, e.g. [18] and references therein). As indication of the size and sign of the asymmetry it has become 
conventional to use the value of its second moment at the reference scale Q 2 = 20 GeV 2 , we obtain: 



which is of the right sign and size as to explain the "anomaly" and furthermore has, as expected, a relatively 
small error due to the dynamical assumptions. Previous determinations [5, 15, 16] generally yield a larger 
value of about 0.0010 to 0.0020 and a typical uncertainty of about 100% or even larger [17]. 

In conclusion, although in our global QCD fits [1, 2] no data with especial sensitivity to the strange 
content of the nucleon have been included, our determination of strange parton distributions, in particular 
Eq. (3), is compatible with particularly sensitive data, e.g. those in [5]. Furthermore, these data induce 
an asymmetry in the strange sea which has been evaluated within our dynamical framework and found, in 
agreement with previous results, to be rather small and positive. This being the case, our original strange- 
symmetric distributions should suffice for most applications. The strangeness asymmetry may, however, 
be relevant for some especially sensitive applications; for these cases our results are available on request 6 . 
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